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The Hamiltonian operator for an unbiased array of Joseph- 
son junctions with gate voltages is constructed when only 
Cooper pair tunnelling and charging effects are taken into ac- 
count. The supercurrent through the system and the pumped 
current induced by changing the gate voltages periodically are 
discussed with an emphasis on the inaccuracies in the Cooper 
pair pumping. Renormalisation of the Hamiltonian operator 
is used in order to reliably parametrise the effects due to inho- 
mogeneity in the array and non-ideal gating sequences. The 
relatively simple model yields an explicit, testable prediction 
based on three experimentally motivated and determinable 
parameters. 



I. INTRODUCTION 

When a potential well propagates adiabatically along 
an electron system that is effectively one-dimensional, it 
carries with it additional electron density, and induces 
dc electric current through the system. Such a pumping 
effect has been studied in small metallic tunnel junctions 
in the Coulomb blockade regime Urn If the propagation 
of the potential well is arranged by phase-shifted gate 
voltages as in Rcfs. 1 and 2 and the potential well car- 
ries a quantised number n of electrons then the induced 
current / is related to the gating frequency / by the fun- 
damental relation / = —nef, where e = 1.602 ■ 10 -19 
C. These Coulomb blockade pumps transporting single, 
normal-state electrons hayerfeached accuracy suitable for 
metrological applicationsM 

Until lately, mainly pumping of single electrons was 
studied, but in a recent articles a quantitative theory of 
pumping Cooper pairs in gated one-dimensional arrays 
of Josephson junctions was presented. It was shown that 
quantum effects render the Cooper pair pump inaccurate 
in case of arrays with small junctions thus explaining 
the failure to demonstrate accurate pumping in the first 
(and so far the only) reported experiment of pumping of 
Cooper pairs B 

In this article, using renormalisation methods, we gen- 
eralise the results derived in Ref. 4 for an unbiased ar- 
ray of Josephson junctions with gate voltages when only 
Cooper pair tunnelling and charging effects are taken into 
account. We consider the supercurrent, the higher order 
corrections for the pumping inaccuracy in case of homo- 
geneous arrays, inhomogeneity of the array, and nonideal 
pumping sequences. Using three experimentally moti- 
vated and determinable parameters we are able to derive 
an expression for the pumping inaccuracy which can be 
directly compared against experimental results. 



It should be stressed, though, that the model is simple 
and neglects many, possibly important features such as 
quasiparticle tunnelling, the coupling to the electromag- 
netic environment and the dissipative effects induced by a 
non-zero bias voltage. On the other hand, without these 
simplifications the problem could not be solved at the 
moment. The efficiency of the renormalisation methods 
in case of Cooper pair pumps is explained by the consid- 
erable symmetries of the model Hamiltonian. RenorraaL 
isation is widely used in atomic and nuclear physics ,ETt£l 
e.g. for the creation of effective two-body interactions 
used in nuclear Shell-Model calculations. 

The article is organised as follows. In Sees, [n] and III 
the renormalisation method is explained and the expres- 
sions for the Hamiltonian as well as the charge trans- 
ported by pumping are derived, respectively. In Sec. [TV 
homogeneous arrays are examined in detail and in Sec.lV 
the inhomogeneity is introduced. Finally in in Sec. VI 
nonideal gating sequences are considered and the theo- 
retical prediction for the pumping inaccuracy is explained 



thus leading to the conclusions in Sec. VII 



II. THE RENORMALISATION METHOD 
A. Few-state dominant systems 

Renormalisation methods may be applied most effec- 
tively in case of "few-state dominant" systems. In such 
systems some of basis states are separated from all the 
others by an energy AE which is large as compared to the 
coupling between any two states in the system. Thus cer- 
tain eigenstates of the full system can be approximated 
by the eigenstates of the " few-state dominant" part only. 
Renormalisation can be described as a bridge spanning 
across the intervening gap. 

For a given orthogonal basis {\l)} the Hamiltonian op- 
erator H can be written as 



H = H + V, H \l) = ei\l), 1 = 1,2,. 



(1) 



where V is the residual interaction defining the coupling 
between eigenstates of Hq . The matrix elements of V can 
often be considerably suppressed by a proper choice of the 
basis states. Explicitly, a part of the Hamiltonian can be 
called " fc-state dominant" if there are k basis states that 
satisfy the conditions 



|ei-em| > A£, 
\V nn .\ < AE, 



if I < k < to, 
for all n.n . 



(2) 
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The requirements of "few-state dominance" are graphi- 
cally depicted in Fig. [l] showing some of the energy levels 
for two examples of a 5-state dominant system. Note 
that Knax, the magnitude of the largest element in V, 
can be much smaller than the energy spread of the low- 
lying basis states. 



B. The effective interaction 

The aim of renormalisation is to generate an effective 
interaction V for a small, active space (P-space) which 
yields the same eigenvalues and eigenstates as the full 
Hamiltonian operator H. Similarly effective operators 
can be defined but this complication can be avoided in 
case of Cooper pair pumps. The renormalisation method 
in the context of nuclear physics is carefully explained in 
Refs. 6 and 7. 

The effective interaction V for a fc-dimensional P-space 
spanned by {|Z)}?=i can be derived as follows. One starts 
from the full Hamiltonian equation 



H\i>) = E\il>), 



(3) 



where E is an eigenvalue, inserts the expansion = 
ai\l) in order to obtain the set of linear equations for 
the coefficients {ai}. In the first k equations of Hamilto- 
nian the rest of the equations can be repeatedly applied 
in the form 



i(>fc) 



E 

i=i 



V rn iai 

E — e m 



Vmm' dm' 



\ ~* * mm' ^m 



V E - £ r 
n'(>k) 



The renormalisation eliminates coefficients ffl m (>fc) 
converges certainly if 



E 

n'(>k) 



| Knm' | 
\E - Cm 



< 1 



(4) 



and 



(5) 



for all states \m) with m > k. Inclusion of additional 
basis states in the P-space improves the convergence by 
reducing the number of coefficients that must be elimi- 
nated. Violation of the property (g) does not necessarily 
imply divergence, but the convergence should be checked 
more carefully. 

We define operators P = Yli=i and Q = 1 — P 
that project onto the P-space and the rest of the full 
space (Q-space), respectively. The renormalised Hamil- 
tonian can then be written as 



H = PH Q P + V. 
where the effective interaction is given by 

Q 



V = P 



n=0 



E - QH Q 



VP. 



(6) 



(7) 



When V converges E is also an eigenvalue of renormalised 
Hamiltonian H which is manifestly consistent with the 
full Hamiltonian equation (||). But, unless V can be eval- 
uated (analytically) up to the infinite order, the renor- 
malisation should be truncated at a given order or when 
preset convergence criteria are met. This approximate 
effective interaction replaces V in Eq. (^). 

The only remaining question is how the eigenenergy E 
should be chosen? The obvious, correct answer is "Insert 
E, get E" . The effective Hamiltonian H yields, allow- 
ing for convergence, k eigenvalues {Ej}^ =1 . By repeat- 
edly inserting an eigenvalue Ej into the definition of the 
effective interaction V we obtain the best available ap- 
proximation for the eigenvalue -Eui of H as well as the 
renormalised eigenstates A good initial guess is the 
eigenvalue E y - ' of the pure P-space Hamiltonian PHP. 
We call this choice " individual" because the renormalisa- 
tion has to be performed for each eigenstate in P-space 
separately. 

The eigenstates {\j}} are not orthogonal in P-space be- 
cause they correspond to different effective Hamiltonians, 
but if they are reexpanded to the full space the orthog- 
onality can be regained. However, the requirements for 
" few-state dominance" ensure that the important advan- 
tage of orthogonality can be used in exchange for only a 
small loss in accuracy. In the " average" choice renormal- 
isation, valid for "few-state dominant" systems, we use 
the average eigenenergy E — k^ 1 Y^=i Hjj m place of 
individual eigenvalues. Thus the evaluation the diagonal 
matrix elements {Vjj} suffices up until the final step of 
the iteration when H has to be diagonalised. 

If the self-consistent iteration is not used we obtain the 
results directly corresponding to initial guesses. We refer 
to these cases as " individual-O" and "average-0" choices. 
These choices have been used since they often offer a 
more transparent interpretation of the result as well as 
much desired analytical results. 

Finally we must emphasise that renormalisation and 
full diagonalisation in a restricted basis are just two sim- 
ilar although unidentical approaches to the eigenvalue 
problem. In renormalisation the full problem is projected 
onto a smaller space while in diagonalisation the prob- 
lem is truncated by discarding all basis states outside the 
restricted basis. 



III. THE COOPER PAIR PUMP 

A. General properties 

An array of Josephson junctions with gate voltages, 
a Cooper pair pump (CPP), in the Coulomb blockade 
regime is an excellent examtjle of a few-state dominant 
system. In an earlier paperu the leading component of 
the inaccuracy for transferred charge in a homogeneous 
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Cooper pair pump was derived. The higher order correc- 
tions to this leading order result are relatively insignifi- 
cant in the immediate pumping regime. In that article 
the inaccuracy was evaluated from the variance of the 
number of Cooper pairs on an island far away from the 
island where most of the charge transfer occurs. 

A crude version of the renormalisation process, 
amounting to "average-0" choice was also used in Ref. 
4 for crosschecking the results. An analytical power ex- 
pansion of the inaccuracy derived below was known as an 
"intelligent guess" but not included in the results. In ad- 
dition, an analytical result for the pumped charge along 
a circular path for N = 3 on the gate voltage plane was 
derived by renormalisation. 

Figure |] shows a schematic drawing of a gated Joseph- 
son array of N junctions. Each junction has a capaci- 
tance Ck and a Josephson energy Ej^. The phase differ- 
ence over the array is <f> £ [0, 2irN). It may be controlled 
by an external bias voltage V over the array according 
to the relation 



(-2e)V 



dt 



(8) 



The oscillation frequency of 4> is approximately V [/iV ] • 
0.5 GHz, but all calculations in this article are done un- 
der an assumption of ideal zero bias (V = 0) yielding a 
constant </>. Each gate voltage V gt k binds charge C g ^kV g ,k 
on island k. In order to (hopefully) transport exactly 
one Cooper pair through the array the gate voltages are 
operated as depicted in Fig. ^. Thus the cycle consists of 
N legs and during k th leg we expect one Cooper pair to 
tunnel through junction k. 

There are two important energy scales in a CPP. The 
first one is the typical Josephson coupling energy Ej re- 
lated to the Cooper pair tunnelling through the junctions 
and the second one is the charging energy Eq related 
to the charging effects of the small islands between the 
junctions. Both are explicitly defined below. The most 
important parameter of the model is the ratio Ej/Eq 
which will be denoted by ej. 



B. The Hamiltonian and supercurrent operators 

When we neglect the quasiparticle tunnelling and other 
degrees of freedom the model Hamiltonian is given by 



H = He 



Hi 



(9) 



where He is the charging Hamiltonian and Hj describes 
the Josephson tunnelling of the Cooper pairs. For an 
array of N junctions the tunnelling Hamiltonian has the 
form 



N 



H, 



cos <p k 



(10) 



fc=i 



where <f>k is the phase difference over the junction k, cor- 
responding to a supercurrent operator 

(-2e)£ J|fc . , -2edH 3 
Ism = i — -sm0 fe = — - — -z-. (11) 



The charging Hamiltonian Hq is diagonal in the ba- 
sis formed by the charge eigenstates \n) where ft = 
{rti, ri2, . . . , n^v-i} and rij is the number of Cooper pairs 
on each island of the array. The normalised gate voltages 
q = {q 1 ,q 2 ,... 1 q N -i} where q k = V g . k C e . k / (-2e) may 
be considered as parameters in He- The diagonal matrix 
elements are given by the classical charging energy 



Ek = E. 



{n,q) 
ch 



JV 

E 

fc=i 



2Ck 



(12) 



where Qk = {—2e)vk is the charge across the junction 
k. We define the typical charging energy of the array 
as Eq = (2e) 2 /2C where the "average" capacitance C 
is given by C = N/ £f =1 C" 1 so that C k = c k C. The 
condition for ideal biasing yields X^fcLi v k/c k — and the 
conservation of charge on each island requires that 



Vk - v k +i = u k 



(13) 



where u = ft — q. The unique solution satisfying these 
conditions is given by Vk = Vk + U where 



N-l 



N-l 



N _ 



fc=i 



j=k ]=1 

By substituting the solution ( jl4| ) into Eq. (|l^) we find 

v 2 



Efl = Ec 



N 2 i / N 

V^- — IV — 



k=l 



Ck 



\k=l 



Ck 



(15) 



where Vk = Vk +y for arbitrary y since the expression ([la), 
is invariant under the transformation {v k } — > {vk + y}El 
This symmetry can be effectively applied when renor- 
malised matrix elements are evaluated. 

The optimal basis for calculations is {|n, </>)}, the basis 
of charge eigenstates augmented by the total phase differ- 
ence over the array <fi = X)fcLi 4>k which is is periodic over 
2Nn. The completeness of the basis was shown in case pf, 
normal-state electron systems by Ingold and NazarovEj 
who also give the canonical transformation between vari- 
ables describing TV separate junctions and variables de- 
scribing N — 1 islands and the array as a whole. The 
conjugate phases on islands {9k} as well as the average 
number of tunnelled Cooper pairs Af = (1/JV) 53fc=i TO fc 
(rrik naturally corresponds to junction k) are completely 
undefined for the chosen basis states. 

From now on we consider only the case when the phase 
difference over the array cf> is kept fixed by ideal biasing 
and since ^> is a constant of motion for the model Hamil- 
tonian (^J) we can explicitly write 
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JV 



H = H c (q) 



H.c). (16) 



n.k—l 



Here the tunnelling vector 6k describes the change of n 
due to tunnelling of one Cooper pair through the k th 
junction. The non-zero components of 5k are (if applica- 
ble) (5k)k = 1 and {5k)k-i = — 1- Each tunnelling in the 
'forward' direction is thus associated with a phase factor 
e «0/AT rp ne corresponding supercurrent operator is given 
by 



Is,k — 



{-2e)E, 



I. A: 



2h 



EHI™ + 4)<n|e^ JV + H.c.) (17) 



where I Ct k = (—2e)Ej t k/h is the critical current of junc- 
tion k. We also define the (average) supercurrent opera- 
tor I s by 



1 N 



(-2e) 8H 



(18) 



The second equality follows from the ^-independence of 
Hq and the relation between (f>k and <f>. The common ex- 
pectation value of Is and 7s, fc hi a stationary state |m) is 
given by (— 2e/h)dE m / d(j) where E m is the correspond- 
ing eigenenergy. Since Is can be expressed as a derivative 
of the full Hamiltonian operator its matrix element be- 
tween two different stationary states |m) and |Z) can be 
expressed simply as 



{m\h\l)<i> 



{-2e){E m -E l ) 



lim 



HZ), 



(19) 



This expression is well-defined although one must keep 
track of the physically unimportant total phases of the 
wave functions. 



C. The supercurrent and the transferred charge 

As shown in the previous articled there are two mech- 
anisms of Cooper pair transfer in the array. The first 
one is the direct supercurrent flowing through the whole 
array due to non-zero <f> and the other one is pumping, 
the charge transfer in response to the adiabatic variation 
of the injected charges q. 

The expressions for the charge transferred by these 
mechanisms can be derived as follows. For each instant 
of time t we introduce the basis of instantaneous eigen- 
states { mrf))} with eigenenergies {E m } of the full Hamil- 
tonian (|) for a given q(t). Assuming slowly varying gate 
voltages we may solve the time-dependent Schrodinger 
equation with the initial condition \ip(to)) = |W(j )) to 
obtain 



IV>(i 



+St) 



_ -iE m St/h\ 



Hto) 



,-iEiSt/h _ e - iE ^ St / h ^Ei\V ^ - ^ 
\ m (t )) + \ Sm (5t) 



E 



t(E l - E m )/h 



at I/ 



(20) 



Here the term | Vj) • §| is the directional derivative of the 
ground state with respect to the change in gate charges q. 
The amount of charge that passes through the junction 
k during a short time interval 5t is then 



SQ k 



i-to+St 

III! 

f2Re 



{lp{t)\hA^{t))dt = St(Is,k)\ m(t0) ) 
■to+St 

(m {to) \Is, k \5m (t _ to) )dt 



In 



(21) 



where we have neglected the term quadratic in \Sm) and 
oscillatory terms by assuming that St ^> h/(Ei — E m ) 
holds for all I. The first term gives the charge transferred 
via direct supercurrent. The second term, the induced 
charge transfer, can be integrated yielding 



-2K Im 



(m\I s , k \l)(l\6m) 
Ei — E„, 



(22) 



where \6m) is the change in the instantaneous eigenstate 
induced by the change q(to) — > q(to + 5t). 

For a closed path 7 the transferred charge must be 
equal for all N junctions so it can be written in terms of 
the average supercurrent operator Is- The total amount 
of charge, Q, transferred over a pumping period r is given 
by 



Q_ 

-2e 



1 



dE m {t) 



dt 



7 l&m) 



Im 



(m\l s \l)(l\dm) 



Ei - E„ 



(23) 



where \dm) is the differential change of \m) due to a dif- 
ferential change of the gate voltages dq. It should be 
noted that the pumped charge depends only on the cho- 
sen path while the amount of charge transferred by direct 
supercurrent also depends on how the gate voltages are 
operated on the path. According to Eq. ( |l9| ) the pumped 
charge for the state \m) simplifies to 



Q P 

-2e 



= 24,^ Im 

7 l(jtm) 



lim 



-(l\dr 



(24) 



Thus the pumped current is mediated by the induced 
mixing of other components into the initial state and 
modified depending on how the relative phases of the 
eigenstates change with respect to the diffential change 
in <j). This formulation proves to be especially effective 
in the regime of two-state dominance. The operating 
frequency of gate voltages must satisfy / Ej/h so 
that the adiabatic approximation is valid, though. 
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D. Numerical, renormalised and analytical results 

In the following sections we refer to our results as nu- 
merical, renormalised or analytical. Numerical results 
are obtained by diagonalising the Hamiltonian operator 
(0) in a given basis and using the corresponding eigen- 
states in order to evaluate the required observable. The 
pumped charge Q p is obtained by numerically integrating 
the second term in Eq. (p3|). 

Renormalised results are obtained by a semianalytical 
process where the renormalised matrix elements in Eq. 
(0) for a given E are expressed analytically but the it- 
eration process is naturally done numerically. Although 
restricting the analytical renormalisation into a given ba- 
sis is quite difficult it facilitates direct comparison be- 
tween the renormalised and numerically obtained results. 
Purely numerical renormalisation forfeits so much infor- 
mation that it is not used in this article. 

The analytical results are obtained by renormalisation 
in such a manner that they can be expressed in a closed 
form. The analytical results for the supercurrent are 
given by the relation (—2e/h)dE m /d<j> while the pumped 
charge is evaluated using Eq. (p4|). 



IV. THE HOMOGENEOUS COOPER PAIR PUMP 

A. The properties of the charging Hamiltonian 

For uniform arrays all Josephson energies are equal to 
Ej and Ck = 1 so the charging energy is given by 



Ex — Er 



N 



provided that Vk are solutions of Eqs. (JT^) . We shall now 
examine the properties and the symmetries of the charg- 
ing Hamiltonian Hc{q). Since the number of Cooper 
pairs on each island may obtain only integer values the 
degree of symmetry of H c (q) depends on how exactly the 
gate voltages q match the set {n}. 

The most symmetrical case is obtained when q = Hq 
for some charge eigenstate |no)- We then find Efi = 
and 



E no±{g) =E c s(N- 



s)/N 



(26) 



where (s) consists of s jumps through s different 
junctions in "forward" direction. Thus the Cooper 
pair pump is single-state dominant near the point 
of high symmetry and the supercurrent is given by 
(— 2e/h)dVii/dcf>. The "average-0" choice now yields a 
supercurrent Ic cos(4>)(Nej/2) N ~ 1 N/(N — 1)! in agree- 
ment with numerical results at q — uq. 

The degree of symmetry is almost as high on any line 
uq + x ■ S r for any r and x £ [0, 1]. All junctions except 
the junction r are equivalent with respect to the charging 



energy. In the region x ~ 0.5 the system is manifestly 
two-state dominant. This limit is the most relevant one 
for us since it is realised by ideal saw-tooth gating shown 
in Fig. |. 

The Hamiltonian is symmetric also at the so-called res- 
onance point q = no + . . . , 1/N), where the states 

{Ko + Sj=i^)}feLi become degenerate. For N degen- 
erate levels with nearest-neighbour coupling the ground 
state supercurrent is given by lies (0) = Ic sm((j)/N)/N 
with <p G (— 7r, 7r) exhibiting a cusp at <p = ±7r. Because 
£,i > the coupling to the other states enhances t he su- 
percurrent which will be explicitly evaluated in Sec. IV C . 



B. Pumping and supercurrent for homogeneous 
arrays 

We will now evaluate the inaccuracy in the pumping 
for the uniform array when the gate voltages are oper- 
ated as depicted in Fig. ||. Due to the symmetry of the 
charging Hamiltonian Hq it is enough to consider any 
one of the legs and multiply the results by N. In the 
Coulomb blockade regime for the saw-tooth gating cy- 
cle the system is always dominated by either one or two 
charge eigenstates. The pumping mainly occurs when 
these two states for each leg are nearly degenerate. 

A two-level Hamiltonian can always be decomposed as 



H = 



ve -i8W 
£2 



(27) 



For the truncated two-level system we have 9{4>) — <f>/N, 
Ej/2 and tj = E^ , j = 1,2. The proper decom- 



(25) position of the renormalised Hamiltonian is 



H 



E, 



(i) 

ch 



iM-.E c v\b{4>)\e-^/ N+ ^ 



v\b(()>)\e i W N+ <t">) 



E 



(2) 
ch 



where v = —Ej/2, a 



+ a 



,(2) 



0'). 



E C 



(28) 



1,2, 



b = bo + b-\e~ 



i<t> 



bie 1 ^ and e 



i<t>b 



b/\b\. The leading 



components for these coefficients are ao oc ej, a\ oc £j 



N 



b 



1 



ce 



i. N- 
0-1 OC £j 



and ti«£j. Thus H clearly 



tends to the unrenormalised Hamiltonian H in the nat- 
ural limit £j — * 0. 

The actual values for these parameters are discussed 
below. The next corrections in H are 02 and b^2 which 
are further suppressed by a factor . The decompo- 
sition ( p8| ) is valid also for inhomogeneous Cooper pair 
pumps with only superficial changes. The renormalisa- 
tion process itself is more complicated, though. 

In the context of this two-level model the renormal- 
isation coefficients a and b have the following interpre- 
tation. The diagonal coefficient ao corresponds to those 
tunnelling sequences that end in the same state inside 
the active P-space and without transporting any Cooper 
pairs through the array. Coefficients a\ correspond to 
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those sequences that transport one Cooper pair in for- 
ward or backward direction, thus yielding the cos In- 
dependent term. The coefficients b arise from the se- 
quences that connect the two charge eigenstates via the 
Q-space, transporting —I, or 1 Cooper pairs through 
the array. Each intermediate state \n) naturally intro- 
duces an energy denominator E — Ej% . 

This interpretation is illuminating and extremely help- 
ful when evaluating the renormalisation coefficients but 
it has a very severe drawback. The picture we obtain 
is, unfortunately, false. The quantum mechanics implies 
that if e,j ^ all charge eigenstates simultaneously co- 
exist although the amplitude of most of these states is 
negligible in the low-lying eigenstates of the system. The 
renormalising sequences appear when we take into ac- 
coun t th e existence of high-lying states as described in 



Sec. [IB 



In order to calculate the integral (|24|) for the leg in 
pumping we will use a parameter rj = (ei — £2)/2v which 
is linear in the ascending gate voltage for the truncated 
system and almost linear for the renormalised system. 
The gating induced correction for the ground state yields 
a term (2|dl) = \dr\j (I + r; 2 ) which is real. Thus we only 
need the imaginary part of the limit in (E4h which reads 



Im 



lim 



0o( 1 l 2 >0o+d0 



d8/d(f) 



2 VI + V 2 



(29) 



For the truncated system dd/dcj) = l/N and the pumped 
charge is given by 



-2e 



1 

2N 



(30) 



The symmetry of the full Hamiltonian (g) implies that 
the pumped charge for the full cycle should be exactly 
— 2e so we bluntly assume that the charge transfer in the 
limit 9 — > (j>/N is exactly Q p = — 2e. We can partially 
justify this assumption by allowing for the missing charge 
transfer via higher excited states and noting that the 
identification of the initial and final states changes after 
each leg. 

For the renormalised system dd/d<j) may be evaluated 
analytically yielding 



d9 rcn 1 



bo(—b-i + bi) cos(f) - b 2 



N 



H<PW 



(31) 



The pumping inaccuracy is then given by a weighted av- 
erage of Nd0/d<p on a single leg. The weights can be 
obtained from Eq. (|3(])but for practical purposes it suf- 
fices to evaluate Eq. ( J3f| ) at the degeneracy point. The 
coefficients are obtained by using the "average" choice 
for the eigenenergy E. In most cases the renormalisation 
includes all terms up to the third order and coefficients ai 
and bzpi up to and including order . The leading cor- 
rection from 6_icos0 is proportional to e^~ 2 as shown 
in Ref. 4. 



In Fig. |j| the pumped charge Q p for N — 3 is studied 
as function of the phase difference (f>. The renormalised 
and numerical results are in good agreement and they 
clearly indicate that the deviations from the leading order 
result (Q P /(— 2e) = I — 9£jcos</>) are important even 
when £j is relatively small. The maximum value for the 
numerical and renormalised pumped charges are 2.14 and 
2.18, 2.87 and 2.98, 3.70 and 3.89 for e } = 0.1, 0.15 and 
0.2, respectively. 

Both minimum and maximum values of Q p correspond 
to vanishing supercurrent which suggests that, in princi- 
ple, the phase differences (f> = and 4> = tt can be dif- 



ferentiated and the ratio Q 



p. max /Q P .linn 



could be used in 



order to determine e j . More realistic models are required 
in order find out if this signature can persist when effects 
due to the electromagnetic environment are included. 
In units Eq the renormalised eigenenergies read 



Ei 

E 2 



<u 



-T 1 - 

' 2 



(Ae)2+ e 2|&(»f 



(32) 



where Ae = £1 — 62- The ground state supercurrent is 
obtained by deriving E\ with respect to <j) with result 



(^s)g.s. = Cosine 



(2) 



2ej 



bo(b-\ + 61) + 26_i6i cos0— (a\ 



(i) 



,(2) 



)Ae/ej 



(2/ £J V(A6) 2 +^|&(^ 



(33) 



In Fig. [| the renormalised and numerical supercurrents 
are shown for N — 5. The renormalisation using the "in- 
dividual" choice reproduces the supercurrent well in all 
three cases. For clarity a basis with 40 states was used 
although it is not large enough to produce the leading 
order (e^ -1 ) supercurrent fully. The slight underesti- 
mation of the supercurrent at the degeneracy point is 
explainable since we could use only leading order terms 
in the renormalisation. Below we study various effects 
related to restricted bases. 



C. The basis-dependent effects 

In order to reliably evaluate the pumped charge or the 
supercurrent one must first select a proper basis in which 
the calculations are performed. The basis should be as 
large as possible so that the discarded states are not im- 
portant but on the other hand such calculations may re- 
quire prohibitive amounts of CPU-time. Our aim is to 
circumvent these problems by using renormalisation tech- 
niques. 

Numerical calculations have been mainly performed in 
three classes of bases we refer to as a-, b- and c-bases. 
From here on we reserve the superscript in parenthesis 
( fe ) for the ^'-dependent part of any coefficient and the 
distinction between charge eigenstates 1 and 2 is taken to 



G 



be complied implicitly. An a-basis contains all states con- 
tributing to the leading order inaccuracy (6^ 2) ) while a 
b-basis produces the leading component of the supercur- 
rent (ai ) fully. Even larger c-basis contains all states 
required for the next-to-leading correction of the inaccu- 
racy (b^ , b[ N \ 

For each length of the array N these bases can be cre- 
ated as follows. The leading component of the inaccu- 
racy (supercurrent) is carried by the N — 1-step (A^-step) 
"paths" containing at most one tunnelling through each 
junction that connect the initial state to the final state 
(itself) for each leg. The total number of necessary states 
is 2 N ~ 2 N and 2 N ~ 1 N for an a-basis and a b-basis, re- 
spectively. A short reasoning confirms that a state should 
be included in a c-basis if it can be reached from some 
state in the corresponding b-basis by a single tunnelling. 
The number of states has not been generally resolved but 
in cases A~ = 4 to N = 8 the c-basis contains 100, 325, 
966, 2695 and 7176 states. 

The numerical evaluation of the pumping integral in 
Eq. ( p3| ) may require up to 10 5 to 10 6 complete diago- 
nalisations per leg before convergence is obtained which 
is extremely time-consuming for bases containing more 
than few hundred basis states. The supercurrent has 
been evaluated even for the A" = 8 c-basis since only one 
eigenstate is needed. Due to computational necessities 
some modifications of the above-mentioned bases have 
been used. 

The differences between bases can be illuminated by 
performing an "average-0" choice renormalisation at the 
degeneracy point of the saw-tooth gating path. Inserting 
the coefficients bj (j = 0, ±1) and respective powers of 
£j into Eq. (|3l"| ) one obtains a power expansion of the 
inaccuracy for small values of ej. This expansion has to 
be corrected for the drop in ground state energy induced 
by the terms Oq . 

In order to include all contributions up to the next-to- 
lcading order we need the expressions for coefficients 



The power expansion of the pumped charge then reads 



(2) >,(2) h (N 2) ^(iV) frj . Simple expressions are 



b w b { 
°o ' -1 



obtained for 



(2) N-2 N(N-l) 
"o,a 4 + 4 (7V _ 2 ) ' 
N-l N(N-l) 



i {2) 

0,b+ 
a(2) 



N 



4(AT-2) 4(2JV-2)' 



= N/2, 



6 (2) = 
°0,b+ 

u(N-2) _ 
J -l,a+ — 



N(N- 1)/2(AT- 2), 



N\ N - 2 N 



(N — 2)! 



where index a and b corresponds to a- and b-bases and 
+ implies that coefficient does not change when basis is 
enlarged. The analytical expressions for the coefficients 



b_i and b\" 1 are composed of several multiple summa- 
tions. The obtained values of b^ for different bases and 
for c-basis are given in Table |. 



Qp „ 

-2e ~ 
+s 2 (t 



COS( 



,(N-2) 



(Nl(N)a^+b^)b^ 



(34) 



where l(N) rj 1 stems from the energy denominators. Its 
value is 1, 1, 11/12, 5/6 and 137/180 in cases N = 3 to 
N = 7, respectively. For A" = 3 and N = 4 the strong de- 
viations from cos((/>)-dependence are explained by addi- 
tional terms 27(ej cos 0) 2 — 81(e,j cos 0) 3 and 24£j cos 2 <j>, 
respectively. The expansion ( |34| ) for N — 3 does not 
compare too well against numerical results in Fig. ^| but 
inclusion of the above-mentioned terms improves agree- 
ment considerably up to ej rj 0.1. 

Next we take a closer look at the case A^ = 5 in Fig. |5| 
where the power expansions foipthe b-basis and the 240- 
state basis (almost full c-basisllia) as well as the results 
for "average" renormalisation are compared to numeri- 
cal results for = 0. The inaccuracy is given in units 
£j which allows more detailed comparison of the pre- 
dictions. The renormalised values follow the numerical 
results more closely than the power expansions but the 
differences between bases are still reproduced well up to 
ej Ri 0.1. In addition the inaccuracy for a-basis is cor- 
rectly placed in between these bases. 

Similar overestimation can be seen in the inset of Fig. |5| 
showing the numerical and renormalised inaccuracies for 
N = 7. Although the results may not seem to be so good 
at the first glance, one should bear in mind that the 336- 
states basis for which the convergence is the best corre- 
sponds to even smaller an inaccuracy than the N = 7 
b-basis. Actual inaccuracy should be evaluated for much 
larger c-basis which is, unfortunately, clearly impossible. 
The scaling of the inaccuracy by £j certainly exaggerates 
the error, too. In conclusion we may state that the renor- 
malisation seems to be able to reproduce the behaviour 
of the inaccuracy reasonably well for any A^ and £j in the 
Coulomb blockade regime. 

The enhancement of the supercurrent at the reso- 
nance point has been studied but the conclusions remain 
valid also in its vicinity. For small values of £j next-to- 
nearest neighbour coupling yields approximate supercur- 
rent 1 + Nsj cos(0/A) in our units of choice, 7 r ei (</>)■ In 
the more general case have used semianalytical third or- 
der renormalisation with 2N(N — 1) state P-space and 
compared the results to the supercurrent obtained by di- 
agonalisation. 

For <j) Ri 7T and 6 < A^ < 10 and the comparison is 
shown in Table |l| clearly indicating that the differences 
between bases for £j = 0.1 are not significant but for 
£j = 0.2 they are growing. The renormalisation calcu- 
lations indicate that for £j = 0.1 the convergence is fast 
both with respect to the order of renormalisation as well 
as the basis. As conclusion we may state that the en- 
hancement is important for large A^ and £j although it 
will not cancel the overall suppression ~ 1/A 2 of the 
maximal supercurrent. 
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V. INHOMOMOGENEITY IN THE ARRAY 

In this section we will derive the leading contributions 
for the pumping inaccuracy and the supercurrent on the 
saw-tooth gating path for an inhomogeneous array. Gen- 
eral considerations imply that the quantity EjEc is ap- 
proximately constant for all junctions in an array. The 
Josephson energy is inversely proportional to the nor- 
mal state resistance Rt of the junction and Eq is in- 
versely proportional to the capacitance C of the junc- 
tion. Since i?x is inversely proportional and C is directly 
proportional to the area of the junction, the product is 
approximately constant for different junctions in the ar- 
ray. (This argument works for junctions fabricated in 
the same batch; otherwise the constants of proportional- 
ity are different.) 

The model Hamiltonian (^) is uniquely defined by the 
ratio £j and relative capacitances {c/.}£Lj once we set 
Ej.k = Cfe-Bj in the tunnelling Hamiltonian Hj. We de- 
fine the inhomogeneity index of the array 



X 



inh = "T7 




(35) 



where = 1/cfc— 1 in order to study the behaviour of the 
inaccuracy as a function of Aj n h. Although the definition 
is valid for arbitrary X; n h the limits A"j n h < 0.15 and 
\gje\ < 0.5 are reasonable for the current technology at 
capacitances of the order of 1 fF. 

On r th leg of the gating cycle the initial and final states 
are \n r ) and \n r + 5 r )(\n r +i)), respectively. The gate 
voltages are given by q = n r + xS r where x € [0, 1] is the 
normalised ascending gate voltage so that x = | corre- 
sponds to the degeneracy point. In order to obtain the 
leading order contributions for the inaccuracy and super- 
current we must set clq = 0, bo = c r in the renormalised 
Hamiltonian (|28| ) and evaluate coefficients a% and 

All the required states on the r th leg can be chosen 
from the classes \n r ± (s)) for s = 0, 1, . . . , N — 1 where 



(s) defined in Sec. [TV A| does not contain 5 r . Let a denote 
a permutation of the set {1, 2, . . . , 7V}\{r}, cr(s) the set 
of s first elements in a and Uk the fc th element of cr.O 
Each a(s) then defines two states with charging energies 



^ [(N/Cr ~ l/cf.) X* 



E± s , x = E nr±(S) {x, a) 

+ (s + G%) {{N -s)± 2x{\ + g r ) - G%)\ (36) 

where G s a — ^2k = i 9a k - The charging energies for the 
initial and final states are given by Eq^ x and Eq^i-x, re- 
spectively. 

The leading order of pumped charge can be obtained 
for the "average-0" choice at the degeneracy point for 
each leg yielding result 



■^5- = 1- (A - COB. 
-2e 



N 



'EE 



(37) 



exact in the limit £j 
and 



0. Here K 



( EMa. 

\ 2E C 



JV-2 



AE^^{N/E C ){E_ S ± -E 0ti ) 

= (s + G s cr )(N-s-l-g r -G s a ). 



n N 

llfe=i c fe 



(38) 



The analytical result < p7\ ) gives us the theoretical ra- 
tio between inhomogeneous and homogeneous inaccuracy 
which will be denoted by Wi n h and compared to numer- 
ical results. The interpretation is obvious since for small 
values of ej the higher order corrections are not very im- 
portant and even then their behaviour is relatively similar 
to the leading contribution. 

The inaccuracy is invariant under arbitrary permuta- 
tions of the set {c^}. Numerically the pumped charge 
may be evaluated for any of the N junctions but far bet- 
ter numerical convergence is obtained by using the aver- 
age supercurrent operator /§. The total inaccuracy for 
an inhomogeneous array is always larger than the corre- 
sponding homogeneous array and the junctionwise inac- 
curacy for Cfc > 1 (cfc < 1) is smaller (larger) than average 
inaccuracy. 

In Fig. ^| Winh is plotted as function of Aj nn corre- 
sponding to some specific sets of relative capacitances 
{cfe} for N = 4 and N = 5. The numerical results have 
been obtained for b-bases. The agreement between an- 
alytical and numerical results is good showing that the 
effects due to inhomogeneity of the array can be reliably 
treated as a correction factor when relative capacitances 
Ck are given. 

The effects due to inhomogeneity can be parametrised 
by obtaining limits for W; n h as a function of Aj n h. For 
-^inh — Iff | this is achieved by considering the even dis- 
tribution of inhomogeneity (g dd = 9, ffeven = — 9 f° r 

1/2 



/ \ 1/2 

even N and gr odd = g ) , ffcvcn 



=9{m)' r '»- 



odd N) yielding a lower limit and maximally distorted 
distribution (gi = gN, gk(>2) = g/N, N — \JN — 1) cor- 
responding to an upper limit. The upper limit yields a 
simple, analytical result 

W inh (X inh , N) < max[/(A inh , N), f(-X inh , N)} (39) 

where 

[1 - g/N] 5 - 3N 



f(9,N) 



N(l+gN) 



[(1 + gN) 2 



N-l 



-(i-.9/^) 2 E 



N-2 N-2 

S i r S 



n i^n 



s + 7 J -- L s + 7 

fc=l Ls=N-k ' s=k ' 



(40) 



with 7 = Ng/(N — g). The analytical expression for the 
lower limit is obtained by explicitly inserting the even 
distribution in Eq. ( [§71) and using the symmetry in order 
to reduce the number of terms to be calculated. Even 
simpler a lower limit can be obtained by considering the 
asymptotical behaviour of the inhomogeneity. We find 



W inh {X inh ,N) >l + a% nh) X 



mh 



(41) 
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where the A^-dependent constant a 
ated from Eq. ( 
42317/3600, 40267/3151, 
to N — 9, respectively .113 



(inh) 



JV 



can be evalu- 
yielding values 8, 85/9, 1279/20, 
and 13.769 for cases N = 4 
In Fig. ^ we graphically present 



the limits for (Wi nn — 1)/Xf nh as function of AT in h m cases 
N = 4 to N = 7. For X inh = 0.15 we obtain limits 20 %, 
24 %, 28 % and 32 % as the maximal increase in inaccu- 
racy as compared to the homogeneous case for N = 4 to 
N = 7, respectively. 

The leading order renormalised supercurrent may be 
evaluated using the "individual-0" choice as follows. The 
eigenenergies of the truncated Hamiltonian are 



Ei 

E 2 



T 3 



A 2 



(42) 



where Ao |X = -Eo.z — Eo,i- x - Using E\ in the renor- 
malisation now yields the leading order supercurrent 
I r>x = (Is)/ ria A, on leg r for ascending gate voltage x 
as 



A' 



EE 

a 1=1 



(7 e sin^) (E J /2f- 1 llZ =1 c k 



(it 1 



AS 



("si 

where energy differences are given by AE a = E s>x — E\ , 

for s 
-(N 

From Eq. (|33j) we find AE L 
ing energy difference. In Fig 



(r,m) 



(43) 



1, . . . ,N - 1 and AE { J' s) = E s , x - E u s = 
s — 1) corresponding to the last elements of a. 

! " " 1 — E2 — E\ for the remain- 
the analytical prediction 



(|43j) is compared against and numerically evaluated su- 
percurrent for N = 6 b-basis. Each curve corresponds 
to a randomly chosen set {c k } £j € [0.02,0.06] as seen 
from the different widths of the peak. The numerical 
and analytical results practically coincide. 

Finally, by using (43) we can explain the a-basis super- 
current in Fig. |]. For a homogeneous array = and 
the charging energies E± SyX are identical so the summa- 
tion over <j yields a prefactor (N — 1)!. The omission of 
certain states amounts to disallowing some paths and the 
numerator has to be corrected by a factor (/ — 1)/(N — 1) 
which is exactly how the supercurrent for a-basis can be 
evaluated. 



VI. THE PUMPING INACCURACY AND 
NONIDEAL GATING 



In this section we will derive 
tion induced by nonideal gatinj 
fine below. When all gate voltage 
in beginning of the first leg the 
be expressed as <fm = n\ + <f ff 
errors. The maximum values of 
voltages are given by q sw ,k = 1 
gate voltages on k th leg read 



the leading order correc- 
; sequences which we de- 
es are turned off in Fig. |^ 
actual gate voltages can 
defining the N — 1 offset 
the normalised sweeping 
whence the initial 



ysw 



Q(k>r 



= n k +q off +q^ 1 '>A 



k — l 



(44) 



where A& = J2j=i I n general, the sweeping voltages 
can be determined more precisely than the offset volt- 
ages, and for most of our calculations we have used 1 % 
and 2 % precisions for them, respectively. 

As in the case of inhomogeneity we can determine the 
effects due to nonideal gating sequences provided we can 
evaluate the charging energy differences at the degener- 
acy point. On the r th leg we choose the coordinates of 
the degeneracy point as 



°(deg,r) 



= n, 



N 

■ E 



(45) 



thus defining quantities fij subject to condition /i r = 0. 
The degeneracy condition E.a r — E^ r+1 can be solved 
easily yielding 



N 



E^'/ C i) = °- 



(46) 



Thus on each leg one must find where the line connecting 
the initial and final gate voltages crosses the hyperplane 
defined by Eq. ( fht ) . This clearly implies that the correct 
nonideality parameter is 



X, 




leg=l j=l 



1/2 



(47) 



which can easily be evaluated once the offset and sweep- 
ing voltages are given. 

The leading order inaccuracy may be evaluated us- 
ing the charging energies ( ]T5| ) at the degeneracy point 
and inserting the corresponding energy differences (mul- 
tiplied by N/Ec) into Eq. J37|). Dividing the result by 
the inaccuracy for a homogeneous array and ideal gating 
sequence we obtain the ratio IFnon . We are mainly inter- 
ested in the behaviour of W non for very small values of 
X non so we have only evaluated the asymptotical limit 



(non) Y 2 



(48) 



for homogeneous array. Here a^ 011 ^ is given by 40/3, 
1225/108, 41/4, 258181/27000, 6136/675 in cases N = 4 
to N = 9, respectively. 

In order to show that X non really is the correct pa- 
rameter we chose several sets of offset and sweeping volt- 
ages for homogeneous arrays. We then evaluated Wnon 
both analytically and numerically for = 4 a-basis with 
Ej = 0.05 and N — 5 a-basis with ej = 0.04. The results 
are shown in Fig. ^| which also show the asymptotic lim- 
its. Theoretical values lie on the curves as well as most 
of the numerical data points. 

The full inaccuracy may be approximately understood 
in terms of the contribution from the inhomogeneity and 
nonideality. In order to show this we have used non- 
ideal gating sequences with some of the inhomogeneous 







arrays already used in Fig. ||. We have collected the 
results in Table III which includes the parameters X^ 
and X non and numerical ratio W non as compared to the 
full theoretical correction and product of corrections due 
to pure inhomogeneity and pure nonideality. The reason- 
able agreement between theoretical and numerical results 
shows that we really can take into account both inhomo- 
geneity and nonideal gating sequences. 

Finally we must note that since the connection be- 
tween X non and offset and sweeping voltages is much 
more complicated than the connection between relative 
capacitances and X m h, there is no straighforward way to 
obtain A„ on from the experimental data. An approxi- 
mate upper limit can be given easily, though. A sort of 
worst-case scenario for |g ff,k + 9sw | < x nou when the pre- 
cision of the gating for each component is known, yields 



Xnon (AO 



Znon(14 + 11A + AN 2 + A 3 ) 1 / 2 

6 1 / 2 A 



(49) 



which can be used in Eq. fltq) and combining this result 
with the estimated limits for Xj nn one obtains reasonable 
limits for the ratio Wnon- Multiplying Wnon by the in- 
accuracy corresponding to homogeneous array and ideal 
gating, allowing for indeterminacy of £j, finally yields 
the final prediction of the present model. The predic- 
tion, based on these three parameters, is a range inside 
which the inaccuracy is expected to lie, but it remains 
to be seen whether the electromagnetic environment or 
other effects strongly modify the present results. 



VII. CONCLUSIONS 

We have studied pumping of Cooper pairs for an unbi- 
ased array of Josephson junctions in an environment with 
vanishing impedance. The present model, which includes 
only charging effects and Cooper pair tunnelling, can be 
reliably solved yielding relatively simple predictions for 
the direct supercurrent and the accuracy of the pumping 
of Cooper pairs. 

We have successfully evaluated higher order corrections 
for the supercurrent as well as the pumping inaccuracy 
for ideal gating sequence and homogeneous arrays. The 
effects due to inhomogeneous arrays or nonideal gating 
sequences can be quantitatively treated by defining pa- 
rameters Ai nn and X non and respective correction factors 

W in h and Wnon- 

The parameters ej and A; n h can be experimentally 
measured and the precision of the gate voltages yields 
limits for A non so the present model can give an explicit 
prediction for the expected range of the experimental in- 
accuracy. The theoretical predictions have been verified 
by numerical calculations, but whether the model is real- 
istic enough to give quantitatively, or least qualitatively 
correct results will be ultimately tested in experiments. 
In any case, further theoretical studies using more real- 
istic and sophisticated models should be performed. 
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TABLE I. The oefficients 



b- and c-basis and 



for a- 

&l for c-basis in cases N = 3 to N = 9 obtained by using 
"average-0" choice renormalisation at the degeneracy point. 
Exact values are given as fractions when the value fits into 
the column. 



N 


b% 


b (N) 

-l,b 


b% 




3 


9/4 


57/8 


69/8 


3/4 


4 


63/2 


436/10 


513/10 


5/2 


5 


83.189 


106.44 


125.54 


5.792 


6 


176.78 


217.07 


261.52 


459/40 


7 


339.51 


405.5 


497.62 


20.834 


8 






894.45 


35.781 


9 






1544.9 


59.135 



TABLE II. The maximal supercurrent (<j> ~ 7r) in units 
Itc1{4>) f° r relatively long arrays and strong coupling. For 
different bases the results were obtained by diagonalisation 
and the renormalised value is for third order renormalisation 
and 2N(N - l)-state P-space. 



N 


£.1 


a-basis 


b-basis 


c-basis 


renorm 


6 


0.1 


1.485 


1.488 


1.491 


1.492 


6 


0.2 


1.881 


1.897 


1.913 


1.914 


7 


0.1 


1.589 


1.591 


1.596 


1.601 


7 


0.2 


2.072 


2.087 


2.113 


2.117 


8 


0.1 


1.689 


1.691 


1.697 


1.708 


8 


0.2 


2.256 


2.271 


2.303 


2.316 


9 


0.1 


1.786 


1.788 




1.815 


9 


0.2 


2.435 


2.448 




2.508 


10 


0.1 


1.881 


1.883 




1.923 


10 


0.2 


2.601 


2.621 




2.690 



TABLE III. The ratios W n on corresponding to nonideal 
gating in an inhomogeneous array. The numerical values 
Wnon.num have obtained by numerical integration for TV = 5 
a-basis with £.j = 0.04 or b-basis with ej = 0.03. The values of 
Wnon.rcn and Wnon, prod are obtained by renormalisation when 
inhomogeneity and nonideal gating sequences are treated si- 
multaneously and separately, respectively. 



Ainh 




^non,num 


^non,ren 


W'noiijprod 


0.0131 


0.0377 


1.0178 


1.018 


1.018 


0.0263 


0.0130 
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FIG. 1. Schematic view of two possible 5-state dominant 
systems. On the right-hand-side all low-lying levels are closely 
packed in energy while on the left-hand-side the spread is 
large as compared to Vmax- In both cases the requirements 
for few-state dominance are well satisfied. 
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FIG. 2. (a) A schematic drawing of a gated Josephson ar- 
ray of N junctions. In pumping Cooper pairs, gate volt- 
ages V gi k are operated cyclically. Ck are the capacitances 
of the junctions and C gt k are the gate capacitances, b) A 
train of gate voltages to carry a charge in a pump. Here 
qk = — Cg ! fcV r gl ) ( ;/2e. The dominant state at the turning points 
of gate voltages are also shown. 
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FIG. 4. The maximal value for the supercurrent in units 
7c£j for N — 5. Curves denote renormalised values and sym- 
bols numerical values corresponding to bases with 40 and 325 
states. In case ej = 0.01 the differences between bases can 
hardly be seen even at the degeneracy point. 
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FIG. 3. The pumped charge Q p /(— 2e) as a function of <f> 
for some values of ej and N = 3. Curves denote renormalised 
values and symbols numerical values which were obtained for 
a 41-state basis. Pumped charge is symmetric in cj> and its 
period is 2n. 
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FIG. 5. The inaccuracy of the pumped charge Q p /(— 2e) as 
a function of ej for different bases and N — 5. Curves denote 
analytical power expansions and symbols numerical or renor- 
malised values. The inaccuracy is given in units £j and the 
phase difference used is <f> — 0. Inset shows the corresponding 
results for the N = 7, 336-state basis. 
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Inh 

FIG. 6. The ratios W- tn h between the inhomogeneous and 
the homogeneous inaccuracies from analytical expression (^) 
and numerical calculations as functions of the inhomogeneity 
index Xi n h which is defined in Eq. (|3f]). Numerical results 
were obtained for b-bases with ej = 0.02 and ej = 0.03 for 
TV = 4 and N — 5, respectively. 
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FIG. 7. The limits for the ratio W ln h as a function of JVi n h 
for array lengths TV = 4 to TV = 7. For small values of X[ n h 
Winh ~ 1 + o^r"^ ■ X? nh which allows several cases to be pre- 
sented simultaneously. 
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Normalised ascending gate voltage 

FIG. 8. A three-dimensional plot of the maximal super- 
current for TV = 6 c-basis in units I c ej for several sets {c^} 
corresponding to different X[ n h. The gate voltages are cho- 
sen from the first leg of the saw-tooth gating path. Junction 
capacitances have been chosen randomly as well as the ratios 
£j which lie between 0.02 and 0.06. Solid curves denote an- 
alytical values and discrete symbols numerical values. The 
modifications of the supercurrent are well reproduced even 
for larger inhomogeneities. The dash-dot curve represents the 
homogeneous supercurrent the largest X inb .. 
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FIG. 9. Wnon for homogeneous arrays and nonideal gating 
sequences as function of X n0 n. The renormalised values are 
almost identical to the asymptotical expansions ( |48| ) shown as 
lines. Numerical values agree reasonably well with theoretical 
results. We used a-bases with ej = 0.05 and ej = 0.04 for 
TV = 4 and TV = 5, respectively. 
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